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Abstract 

Recent experimental investigations and physical modeling studies have indicated that 
turbulence behaviors within a liquid jet have considerable effects on the atomization 
process. For certain flow regimes, it has been observed that the liquid jet surface is highly 
turbulent. This turbulence characteristic plays a key role on the breakup of the liquid jet 
near to the injector exit. Other experiments also showed that the breakup length of the 
liquid core is sharply shortened as the liquid jet is changed from the laminar to the 
turbulent flow conditions. In the numerical and physical modeling arena, most of 
commonly used atomization models do not include the turbulence effect. Limited 
attempts have been made in modeling the turbulence phenomena on the liquid jet 
disintegration. The subject correlation and models treat the turbulence either as an only 
source or a primary driver in the breakup process. 

This study aims to model the turbulence effect in the atomization process of a cylindrical 
liquid jet. In the course of this study, two widely used models, Reitz’s primary 
atomization (blob) [1] and Taylor- Analogy-Break (TAB) [2] secondary droplet breakup 
by O’Rourke et al. are examined. Additional terms are derived and implemented 
appropriately into these two models to account for the turbulence effect on the 
atomization process. Since this enhancement effort is based on a framework of the two 
existing atomization models, it is appropriate to denote the two present models as “T- 
blob” and “T-TAB” for the primary and secondary atomization predictions, respectively. 
In the primary breakup model, the level of the turbulence effect on the liquid breakup 
depends on the characteristic time scales and the initial flow conditions. This treatment 
offers a balance of contributions of individual physical phenomena on the liquid breakup 
process. For the secondary breakup, an addition turbulence force acted on parent drops is 
modeled and integrated into the TAB governing equation. The drop size formed from this 
breakup regime is estimated based on the energy balance before and after the breakup 
occurrence. The turbulence energy is also considered in this process. 
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I. T-BLOB MODEL 


The original formulation of the primary breakup “blob” model is based on the stability 
analysis performed by Reitz [1]. The main contribution of this research effort is to 
incorporate the turbulence effect in modeling of the liquid jet breakup. Terms and 
parameters associated with the turbulence behavior shall be derived and implemented 
appropriately to the existing blob model. Hence, it is proper to denote the present model 
as “T-blob” model. In order to account for the both effects of the surface wave 
perturbation and the turbulence motion, the phenomenological approach is chosen in this 
model development. The resulting model should adequately include the combination of 
these two phenomena in the breakup model when the both co-exist in the considered flow 
conditions. However, the model is also capable to primarily describe a single event only 
if the other is neither dominant nor present in certain flow conditions. In such a case, the 
original formulation should be retained. For all aspects, the model should reasonably 
maintain a smooth transient from the fully turbulent to non-turbulent flow regimes. Also 
in this phenomenological method, the time scales associated to the surface wave and 
turbulence characteristics are used to weight these individual effects. The consideration 
of the time scale for this weighting is based on a reason of which the phenomena of the 
smaller time scale would have a stronger influence in the liquid jet breakup process. In 
the course of this study, the framework of Huh et al. [3] in turbulence characterization is 
employed. The length scale and time scale associated with the primary breakup are 
comprised of the ones described in the blob model that represents for the surface wave 
instability and of the turbulence behavior shown in Huh et al.’s model [3]. 

It is proposed that the reciprocal of the primary atomization time scale is equal to the sum 
of the reciprocals of the time scales associated with the wave motion and the turbulence 
characteristics with the inclusion of the respective weighting factors as follows: 



c w and c, are the weighting coefficients with respect to the wave perturbation time scale 
and the turbulence time scale. T t andx p are the turbulence time scale and the primary 

atomization time scale [1]. In the work of Huh et al. the parent drops carry homogeneous 
isotropic turbulence starting at the injection nozzle exit. By further assuming that no 
additional internal turbulence is generated, they were able to obtain an analytical solution 
of the turbulence scale through the use of the well-known k - e turbulence model. So the 
turbulence time scale can be expressed as a function of time: 

T t =t 0 +0.0828t 1.2 

t is the time referenced from the time at which the parent drop is issued at the injection 
nozzle exit. T° is the initial turbulence time scale evaluated from the initial turbulence 

kinetic energy, k° , and its corresponding dissipation rate, e” , at the injector exit. 
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U is the injection velocity of the liquid issued from the injection nozzle, which has the 
length, L, and the diameter, D. Cd, Kc, and s are the discharge coefficient, the loss 
coefficient due to the nozzle entrance sharpness, and the downstream-to-upstream 
contraction area ratio of the injection nozzle, respectively. Equation 1.1 ensures that the 
process of the smaller time scale would have a more significant effect on the atomization 
process. 


The results of the experiments conducted by Faeth at al. [3-4] showed that the drops 
produced from the on-set breakup regime have their sizes in the same order-of-magnitude 
of the turbulence length scale. Furthermore, in his blob model, Reitz set the radius of the 
product drop to be proportional to the wave length corresponding to its fastest growth 
rate. These indications suggest that it is reasonable to choose the radius of the product 
drops as the characteristic length scales for both the wave perturbation and turbulence 
phenomena. It is proposed that the characteristic length scale is expressed as: 

1 c c, r r 

- — =>r=- , 1.4 

r P r w r t (l-c t )r t +c,r w 

r p is the radius of the product drop. The radius r w associated with the wave motion can be 

determined from Reitz’s model. For the value of r t , the formulation of estimating the 

product drop size proposed by Huh et al. is employed. The drop size can be represented 
by a probability density function (PDF), which is given to be proportional to the ratio of 
the turbulence energy spectrum and the atomization time scale. The notion for this 
representation is that an eddy motion with larger turbulence kinetic energy and a shorter 
atomization time most likely causes the drop breakup more frequently than others 
containing the lower energy level and the longer atomization time. It should be noted that 
the authors assume the atomization time scales to be a linear sum of the turbulence and 
wave growth time scales. This time scale, however, is set equivalent to the primary 
atomization time scale defined in equation 1.1 for the present model. By assigning the 
wave number of the turbulence energy spectrum as the inverse of the product drop 
diameter, the PDF for the product size distribution can be shown in a form as: 
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T p is evaluated by equation 1.1. Lt is the turbulence characteristic length scale and is 
modeled using the analytical solution of the k-e turbulence model [7]: 


L t =L° 
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The initial turbulence length scale L° is obtained from the turbulence kinetic energy and 
its dissipation rate as: 
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The constant C in equation 1.5 is determined by the normalization condition imposed on 
the PDF of the drop size distribution as: 

^(■>,=1 1-8 


In the T-blob model, the rate of change in the parent drop radius is taken the same form 
as the blob model: 

— = — when r < a 1.9 

dt T p 

It should be noted that when a < r p the parent drop would no longer strips its mass to 

create the product drop. An example for parent drop diameter change vs. time is given in 
figure 1. 


Parent Drop Diameter 



Figure 1. Parent drop diameter vs. time 
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The velocity and properties of the product drops in the present model are formulated in 
the same manner as in Reitz’s model. In each breakup event, the product drops were 
given the same fluid properties and physical location as the parent. For the velocity, the 
product drops would carry the same velocity, V, of the parent along with extra two 
normal velocity components, v and w, in reference to the trajectory of the parent drop to 
be added, v and w are given by 


v = ]vj tan (0/2) sin 0 
w =|V|tan(0/2)cos<|> 


1.1 


Where 


tan (0/2) = A,AJQ/U 


<j>is chosen at random on the interval (0,27t)and the constant A) is set equal to 0.188 for 
sharp entrance constant diameter tube injection nozzles. 


Along with the inclusion of the turbulence effect on the primary atomization process, the 
subject phenomena will be also considered in the secondary droplet breakup model, 
which will be discussed in the next sections. At any rate, this model requires the initial 
velocity fluctuation quantity of the product drops right after their formation. In order to 
evaluate such a quantity, the energy budget of the parent drop during the primary breakup 
would be examined. This consideration is based on the notion that the overall energy of 
the parent drop combined with the one of the product drops has to be conserved during 
the breakup process [7], One example of the energy budget is shown in figure 2. Detailed 
derivation will be given in the full paper. 



Figure 2. Energy budget during breakup process 
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II. Derivation of turbulence term in T-TAB model: 


In an attempt to account for the effect of the turbulence on the secondary droplet breakup, 
a term, F t , representing a force associated the turbulence behavior is introduced to the 
right-hand side of the original Taylor Analogy Breakup [2] equation: 

m£ = F+F t -kC-dC 2.1 

In order to derive a force term to account for the turbulence behavior in a liquid droplet, a 
force associated with turbulence kinetic motion of an eddy within a liquid droplet is set to 
a function of a turbulence characteristic velocity scale, 0 , and a turbulence characteristic 
length scale, t , of the subject eddy as follows: 

F-*~P,#'A 2.2 

Where 

p, Density of liquid 
A Cross-sectional area of an eddy 

The proportionality in equation 2.2 is implied to account for the effect of ellipticity, non- 
uniform velocities within the eddy and the efficiency of the conversion of the turbulence 
force to the displacement oscillation and distortion of the droplet. Since the subject eddy 
has a turbulence length scale of t . It is logical to define the cross-sectional area. A, to be 
proportional to this length scale as follows: 

a~ e 

Thus, 


Fe** ~ P,<^ 2 2.3 

It is reasonable to assume that the turbulence characteristic velocity and length scale of an 
eddy with the inertial range of the turbulence spectrum are related as follows [5]: 


Where 

v 0 Time-averaged mean of fluctuation velocity 
% Radial integral length scale associated with turbulence 
The fluctuation velocity v 0 can be obtained from the initial turbulence kinetic energy: 
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In a study of the onset of turbulence primary break up, Faeth et al. chose the radial 
integral length scale, X , and the turbulence length scale, C , to be comparable to the initial 
liquid jet diameter and Sauter-mean diameter (SMD) of the breakup droplet. 
Consequently, the force associated with the turbulence process on a droplet is expressed 
as: 
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Faeth at el. have observed that the liquid jet surface fluctuation and the drop size resulted 
from the on-set breakup have the same order of the turbulence length scale. Hence, it is 
logical to also assume that the droplet distortion displacement in the present study has the 
same order of this length scale. This leads to a hypothesis of having the characteristic 
turbulence length scale to be proportional to the displacement distortion x from the drop 
equilibrium position. That is 


Therefore equation 2.7 becomes 


and equation 2.1 becomes 
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Similar to the method of estimating the product drop size in the TAB model, the energy 
balance between the parent and product drops is used to determine the Sauter-mean 
radius of the product drops in the present model. Hence, individual energy forms 
contained in the both parent and product drops are formulated. The turbulence energy is 
also included in this process. Detailed derivation will be presented in the full paper. 


Examples: 

Figure 3 shows the penetration of a diesel jet into a 5MPA chamber and figure 4 shows 
the snap shot of droplet distribution of the atomizing jet at 4.50 msec. More detailed 
comparison for both non-evaporating and evaporating atomizing sprays will be presented 
in the final paper. 
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Figure 3. Penetration length vs. time 



Figure 4. Spray distribution at 4.5 msec. 
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